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Recent calculations predict a possible existence of nuclei with tetrahedral symmetry: more pre- 
cisely, the mean-field hamiltonians of such nuclei are symmetric with respect to double point-group 
Tj 5 . In this paper, we focus on the neutron-rich Zirconium isotopes as an example and present re- 
alistic mean- field calculations which predict tetrahedral ground-state configurations in 108 - 110 Zr and 
low-lying excited states of tetrahedral symmetry in a number of N < 66 isotopes. The motivations 
for focusing on these nuclei together with a discussion of the possible experimental signatures of 
tetrahedral symmetry are also presented. 

PACS numbers: PACS numbers: 21.10.-k, 21.60.-n, 21.60.Fw 



A highly relevant field of investigation of nuclei with 
unusual combinations of the proton and neutron numbers 
is the search for manifestations of exotic symmetries. In 
our previous work we reported Q , that of all the various 
possible non-spherical shape symmetries, tetrahedral and 
octahedral ones may lead to the largest single-particle 
shell gaps. Indeed, these can be comparable to, and 
sometimes even larger than the well-established spheri- 
cal shell gaps. These properties, far from being restricted 
to nuclei within a given range of mass or isospin, should 
in fact be present throughout the nuclear chart. The 
underlying considerations are based on a very general 
analysis of the point groups of symmetries of the nuclear 
mean-fields. They have been confirmed independently 
via the macroscopic-microscopic method of Strutinsky |l| 
and the self-consistent Hartree-Fock techniques . 

The purpose of this Rapid Communication is to sup- 
port the search for tetrahedral nuclei by providing real- 
istic calculations in some selected study cases. 

The neutron-rich Zirconium (Z=40) isotopes provide 
such an appropriate case. Our previous work predicted 
that, among others, nucleon numbers Z = N = 40 and 
Z = N = 70 correspond to large shell gaps for tetrahe- 
dral shapes and thus can be thought of as 'tetrahedral 
magic numbers' 0, Although the yrast line of the 
N = Z = 40 system 80 Zr has been experimentally stud- 
ied above spins of 10H |4j, the detailed information con- 
cerning in particular the electro-magnetic transitions (see 
below) which is required for probing the tetrahedral sym- 
metry has not yet been achieved. On the neutron-rich 
side, Zirconium isotopes up to the N = 70 system 110 Zr 
have been synthesised using the projectile fission tech- 
nique at GSI, The production of such nuclei opens 
up the possibility of studying their properties via isomer 
spectroscopy following projectile fragmentation (e.g. [(|) 
and projectile-fission (e.g. Q) or beta-delayed gamma- 
ray spectroscopy (e.g. To date however, such studies 
have not yet been performed in the very heavy Zirconium 
isotopes and the most neutron rich Z = 40 isotope where 
information on excited states is available is 104 Ztka, which 



was studied as a fragment from spontaneous fission of Cm 
sources in Q. 

In this article we suggest that the heavy Zr isotopes 
are amongst the best candidates to test the hypothe- 
sis of nuclear tetrahedral symmetry: 110 < 112 Zr are in- 
deed predicted by Hartree-Fock calculations to have 
T^-symmetric ground-state configurations, and both 
non-self-consistent approaches and Hartree-Fock calcu- 
lations agree on the prediction of low-lying excited - 
symmetric minima for a chain of isotopes near 110 Zr. As 
outlined above, these nuclei are viable for future experi- 
mental study. 

In the basis of spherical harmonics, the distance of a 
point on the nuclear surface relative to the centre of the 
nucleus can be expressed by 
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where refer to multipole deformations. In Eq. JQ), 
a = {a\^} represents the complete set of deformations 
and c(a) accounts for nuclear volume conservation. Sev- 
eral ensembles of deformation parameters a = {au} can 
generate a geometry characteristic of the group . The 
simplest case occurs when all deformations except 0:32 are 
equal to 0. This simplest form of the tetrahedral defor- 
mation is also likely to be favoured in nature since other 
realizations of the tetrahedral symmetry involve much 
higher order of multipoles, A > 7, cf. Ref. 0. 

The group is characterized by two non-equivalent 
two-dimensional irreps, and one four- dimensional one. 
This is an unique feature in nature, since, apart from 
the spherical symmetry, only the groups and 
(the latter, octahedral symmetry is discussed in Ref. Q) 
produce single-particle degeneracies higher than 2, cf. 
Refs. [Tol Hi) . This high degeneracy pattern, together 
with a relatively low number of irreps favours the ap- 
pearance of large shell gaps. In our previous work |l|, 
we reported that such gaps are predicted in a number of 
different mass regions, and presented various examples of 
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potential tetrahedral minima. 

In the current work, realistic nuclear structure calcu- 
lations are performed within the standard Woods-Saxon 
mean-field approach with the universal parameterization. 
Figure ^ shows the single-particle energies for exotic nu- 
clei with A ~ 110 as functions of the lowest order tetra- 
hedral deformation obtained from this hamiltonian. The 
most striking result in Fig. ^ 1S the presence of a large 
tetrahedral shell gap at Z = 40, comparable to the spher- 
ical Z — 50 gap. Indeed this gap is larger by ~ 3 MeV 
than the spherical one at Z = 40. 

To calculate the total energy of the nucleus, we em- 
ploy the microscopic-macroscopic method whose realiza- 
tion closely follows that of Refs. 0] and 0]. I n this 
approach, the total energy of the nucleus is given by: 
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where E mn£3 _ is the liquid drop energy, here calculated 
as in Ref. |14| ; the latter realization includes in particular 
the nuclear surface-curvature contributions. 

The quantum correction 5E s heiu accounts for the 
nuclear shell structure. It is calculated accord- 
ing to the standard Strutinsky prescription [Hjj. and 
supplemented with the pairing correlation energy, 
5E pa i r = Ep a i r — E pa i r= o, obtained within the pairing- 
selfconsistent Hartree-Fock-Bogolyubov (HFB) formal- 
ism as in Ref. fl6| . The E pa i r term refers to the energy 
of the system in the presence of pairing, calculated by 
using the Particle Number Projection (PNP) technique 
(before variation) , and E pa i r =o is the analogous term ob- 
tained by setting the pairing correlations to zero. The 
average pairing strength constants have been obtained in 
the usual way by fitting to the known nuclear masses of 
nuclei in the studied region. 

The deformation space used includes a priori all possi- 
ble ct\^ up to A = 12, with —A < fx < +A. However, the 
number of multipoles effectively used was significantly re- 
duced since the effect of high- A was found to be negligible 
(by a direct verification). 

The potential energy landsca pes of Fig. |21 top, confirm 
the earlier calculations by fl3l Il7l Irsl Il9j . all of which 
predict that at its ground-state, the 110 Zr nucleus has 
a significant prolate deformation with an axial hexade- 
capole deformation a^o of about 0.10. The calculation 
for the top frame of Fig. [21 has been performed accord- 
ing to the general procedure discussed above, but the 
deformation space included only the quadrupole (/?2>7) 
and axial hexadecapole 0:40 deformation (in the follow- 
ing we employ the standard Hill- Wheeler parameteriza- 
tion for quadrupole deformations). Beside the prolate 
ground-state, we note the presence of a very flat oblate 
pattern together with a small spherical minimum. How- 
ever, the inclusion of the tetrahedral degree of freedom 
has a strong impact on this energy landscape. While the 
prolate minimum remains, the minimum at 02 = is 
lowered in energy by more than 1 MeV and corresponds 
to a tetrahedral configuration with CV32 ~ 0.14. In the 
calculations this minimum appears as pure i.e. all other 
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FIG. 1: Single-particle neutron (top) and proton (bottom) 
energies as a function of the lowest-order tetrahedral deforma- 
tion (= a.32) for nuclei around 110 Zr7o. The Nilsson labels cor- 
respond to the leading basis wave-function in the expansion of 
the single-particle wave function. The eigenenergies belonging 
to the four-dimensional irreps (thus forming 4-fold degenerate 
multiplets) are marked with the full lines and double Nilsson 
labels. The levels belonging to the two non-equivalent two- 
dimensional irreps are marked with dashed lines. 



deformations turn out to be zero as the result of the min- 
imization. 

The often employed technique of calculations using a 
mesh of deformation points (as in Fig. |2J) has the draw- 
back of taking into account only a very partial projection 
(usually less than 4 deformations) out of the full defor- 
mation space. In order to remove this limitation, we have 
adopted a strategy based on a dynamical minimization 
of the total energy: for fixed values of a given defor- 
mation parameter (in our case the elongation along the 
z-axis, fa), the total energy was minimized with respect 
to many other deformation variables using the standard 
Variable Metric Method described in [20|. In our calcu- 
lations the deformation space is composed of the entire 
set of octupole {0:3^; /1 = 0,1,2,3} and hexadecapole 
{a4^; fx = 0, 1, . . . , 4} degrees of freedom, in addition to 
the 7 angle, i.e. 11 deformation parameters in total. As 
mentioned above, the effect of higher order multipoles has 
been investigated and is neglegible in the present context. 

Figure|21shows the minimized-energy cross-sections for 
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FIG. 2: (Color online) Ground-state total energy surfaces 
in 110 Zr as a function of the quadrupole (/3, 7) deformations. 
Top: minimisation over Q40 at each (/3, 7)-point; bottom: sim- 
ilar to the above but minimized over tetrahedral deformation 
0:32- Note that the tetrahedrally-symmetric minimum is lower 
by more than 1 MeV than the spherical minimum. 
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FIG. 3: Total energy as function of quadrupole deformation 
in a chain of Zr isotopes. Negative values of 02 correspond to 
oblate shapes. Each P2 ~ minimum correspond to a tetra- 
hedral configuration. (The pairing correlations were included 
and treated with the help of the PNP technique). 

the ground-states in several isotopes of Zirconium as 
function of the quadrupole deformation /3%. For each 
isotope three minima appear: one at an elongation of 
/?2 ~ 0.35 — 0.40, another one around 02 ~ —0.25, and 
a third one at {5% ~ 0. The latter corresponds to a finite 
tetrahedral deformation of about 032 ~ 0.14. The curves 




plotted in Fig. [3] may provide an estimate of the height 
of the barriers between the tetrahedral minima and the 
other deformed minima. In the case of 110 Zr, this barrier 
amounts to up to ~ 1.5 MeV with respect to the oblate 



minimum. 



TABLE I: HFB energies (in MeV) for various energy minima 
in heavy Zr isotopes, relative to the energy of the tetrahedral 
minimum. Calculations were done with the SLy4 force [2l|. 



Nucleus 



*Zr 



'Zr 



j Zr 



'Zr 



2 Zr 



Tetrahedral +0.00 +0.00 

Spherical +0.22 +0.29 

Oblate -1.57 -1.52 

Prolate -2.07 -1.76 



+0.00 +0.00 +0.00 

+0.39 +0.43 +0.03 

-1.10 +0.07 +0.30 

-0.68 +0.27 +1.01 



Self-consistent Hartree-Fock-Bogliubov calculations 
have also been performed for these nuclei, using the SLy4 
force and the HFODD code (version 2.07) of pf . 
and are reported in Table [I] The same energy profiles 
as shown in Fig. [21 are predicted, with the noticeable ex- 
ception that the tetrahedral minimum in 110 Zr and 112 Zr 
lies lower in energy than the now secondary prolate min- 
imum. Both Fig. |3 and Table [I] show that the influence 
of the shell effects leading to the tetrahedral symmetry 
diminishes when going away from the 'tetrahedral magic 
gap' N — 70. For the lighter isotopes like 104 Zr and 
106 Zr, the tetrahedral minimum lies very high in energy 
and is extremely shallow. 

An important aspect associated with the preceding dis- 
cussion is to find a possibly unambigous experimental sig- 
nature of the discussed symmetry. Here we would like to 
construct a criterion based on the rotational properties of 
a nucleus with the tetrahedral symmetry. A discussion 
of other possible experimental signatures can be found 
in Q. We use a generalisation of the rotor hamiltonian 
often used in molecular physics, cf. e.g. j2^. We define 
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where the first term which is proportional to the square of 
the nuclear angular momentum operator assures that the 
energy vs. spin dependence is approximately quadratic. 
The second, tensor term, is responsible for the tetrahe- 
dral symmetry of the hamiltonian @. The constant 
jfo represents the effective moment of inertia. By con- 
struction, operators T 3i ± 2 are spherical tensors of rank 
A = 3 with the 'magnetic' components /i = ±2; they 
are appropriately symmetrized functions of angular mo- 
mentum operators {I x , I y , I z }. The imaginary coefficient 
^3,2 is a parameter of the model. According to the ro- 
tor formalism such a parameter must be independent of 
operators {I x , I y , I z }, but it may be a scalar function of 
the spin quantum number /. Here we use h^^ = c^^/I 2 , 
with an imaginary constant 03,2- This is a convenient 
parametrisation that allows to keep the third order term 
dominated by the usual, approximately quadratic energy 
dependence on spin. In the following we use C3.2 = i 0.03; 
it turns out that the exact value of this constant has a 
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very minor influence on the electro-magnetic transition 
rates, the latter representing the main concern in this 
part of the discussion. 

The solutions of the Schrodinger equation with hamil- 
tonian for an even-even nucleus are characterized by 
integer spins and thus transform with the help of the 
irreducible representations of the 'simple' Td group (as 
opposed to the 'double' Tj'-group of symmetry that ap- 
plies for the nucleons (spinors) mentioned earlier). The 
group Td has four irreducible representations, three of 
them, denoted CI, C2 and C3 are one dimensional, the 
fourth one, Tl, is three-dimensional. Because of the lat- 
ter, in the strict symmetry limit the rotational spectra 
must contain three- fold degenerate levels (cf. Fig. 0J. 
However, it is well known that quantum systems man- 
ifest zero-point motion oscillations: in the nuclear case 
the most important oscillations of this type are expected 
to be of quadrupole nature. We have estimated, using 
collective model techniques with the Bohr Hamiltonian 
that the expectation value of the quadrupole deforma- 
tion in 110 Zr is ~ 0.05, depending somewhat on the de- 
tails of the model. To simulate the zero-point motion 
effect in the following we therefore perform the calcula- 
tions using a small quadrupole triaxial "contamination" 
term introducing three moments of inertia J x , J y and J z 
which slightly deviate from Jo. These are calculated us- 
ing uniform density cllispoidal shapes with quadrupole 
deformations — 0.05, 7 = 30° and 0:32 = 0.15. 
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FIG. 4: Color online) Tetrahedral rotor energy spectrum; the 
irreducible-representation symbols are marked on the right 
hand-side. Arrows indicate the strongest E2-transitions [in 
units of e 2 (/m) 4 ] depopulating two lowest states of spin 
1=5, as an example. Full lines give stretched, dashed lines 
the non-stretched AI = 1 transitions. For easier legibil- 
ity the degenerate Tl states are marked as split. Observe 
that the non-zero transitions connect only one-dimensional 
irrep states among themselves and the three-dimensional ir- 
rep states among themselves. In addition, there are no E2 
transitions of the type CI — ■> CI, C2 — » C2 nor C3 —> C3. 

The zero-point motion mechanism simulated in this 
way has a very small influence on the degeneracy (de- 



viations of the order of 2 keV) but it has a tremendous 
impact on the electro- magnetic transition probabilities. 
In particular, at the exact symmetry limit the only tran- 
sitions allowed are E2> ones; in the calculations with the 
zero-point motion effect included, the quadrupole tran- 
sitions are stronger by 5-6 orders of magnitude and the 
non-stretched B(E2)s are usually slightly stronger than 
the stretched ones. 

In the adiabatic approximation, the nuclear hamilto- 
nian is a sum of the rotor and the intrinsic energy term 
(cf. Eq. (4-3) of Ref.|23) and it follows that the eigen- 
energies are also sums of two corresponding terms, while 
the transition rates depend on the electric moments cal- 
culated with the help of the intrinsic hamiltonian. The 
latter were calculated here using an uniform charge distri- 
bution with the tetrahedral deformation of 0. 15 and small 
quadrupole deformations cited above. Results in Fig. 0] 
illustrate only the rotor part of the spectrum. Conse- 
quently when reading them one has to bear in mind that 
the energies of e.g. even-spin sequence and of the odd- 
spin sequence may be shifted with respect to one another 
by an amount that needs to be taken from the interpreted 
measured spectrum. The results of this kind could also 
be used as a guidance for a band starting at I — 3 (tran- 
sitions corresponding to I > 3) or for the bands with all 
parities negative. 

Moreover, it should be emphasized, that the reduced 
transition probabilities for (7 — » I + 1) -transitions are 
usually some percentage stronger than those for (I — » 
/ — l)-transitions, so that if e.g. the odd-spin sequence 
is shifted upwards, the inter-band transitions of the type 
5 — > 6, 4 — ► 5, etc. are to be expected. In short: the 
information and prediction content of the diagrams of 
the Fig. 0] type is much richer than what the first-glance 
reading may suggest. 

The partial transition pattern illustrated in Fig. ^dis- 
tinguishes clearly between the decay modes of a tetrahe- 
dral quantum rotor and the one of an ellipsoidal symme- 
try. In the latter case, only the stretched E2-transitions 
are present, connecting exclusively the band-members of 
the same Z^-symmetric rotor irreducible-representations, 
as discussed in Chap. IV of Ref.|24[- But even in the case 
of a strong triaxial minimum competing (7 ~ 30°) the re- 
duced probabilities for the non-stretched E2-transitions 
are nearly two orders of magnitude weaker compared to 
the stretched ones. Moreover, they connect only the 
states of the same D2- symmetry (same common irreps) 
- while in the case of the T^-symmetry, the (Cj — > C)- 
transitions are forbidden, cf. caption to Fig.0| These dif- 
ferences provide an excellent 'yes-no' criterion, provided 
that a sufficient number of experimental transitions have 
been identified. 

In the illustration of Fig. 0] we have selected a con- 
figuration composed of the lowest lying rotational levels 
of the same parity with spins I < 5; a more complete 
analysis will be presented elsewhere. 

In summary: Low-lying tetrahedral configurations 
in Zr isotopes with 64 < N < 72 are obtained with both 
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the HFB and Strutinsky techniques. The only notice- 
able difference between these very different realizations 
of the mean-field comes from the relative position of the 
tetrahedral minimum relative to the ground-state. HFB 
calculations with the SLy4 interaction predict a tetrahe- 
dral ground-state in 110_112 Zr. The barrier heights which 
separate the tetrahedral minima from their prolate coun- 
terparts have also been estimated and we emphasize that 
this region of the nuclear chart is potentially accessible 
for exploration using isomer and /3-delayed 7-ray spec- 
troscopy following projectile fission. 

The spectral properties expected to accompany the 
tetrahedral symmetry in nuclei have been briefly dis- 
cussed in terms of the electro-magnetic transitions in 
a quantum rotor with tetrahedral symmetry. The pre- 
dicted decay pattern is very characteristic: it includes 
a competition between stretched and non-stretched E2- 
transitions and excludes the possibility of connecting 
some well defined groups of levels. This decay pattern 



is clearly distinct from known transition patterns of both 
the axial or nearly-axial rotors as well as strongly tri- 
axial rotors whose spectra contain principally the E2- 
transition cascades. Since in the discussed nuclei the 
tetrahedral-symmetry minima are in competition with 
the ones of quadrupole type deformation, the electro- 
magnetic decay criteria formulated offer a clear-cut dis- 
tinction. 

We believe that the predicted structure of the nuclear 
potential energy and the corresponding minima strongly 
encourages the search for the first experimental evidence 
of existence of these highly exotic nuclear symmetries. 
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